Semi-discrete solitons in arrayed waveguide structures with Kerr nonlinearity 
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We construct families of optical semi-discrete composite solitons (SDCSs), with one or two inde- 
pendent propagation constants, supported by a planar slab waveguide, XPM-coupled to a periodic 
array of stripes. Both structures feature the cubic nonlinearity and support intrinsic modes with 
mutually orthogonal polarizations. We report three species of SDCSs, odd, even, and twisted ones, 
the first type being stable. Transverse motion of phase-tilted solitons, with potential applications 
to beam steering, is considered too. 
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Power exchange among evanescently coupled cores in 
arrays of single-mode nonlinear waveguides has attracted 
a great deal of interest, because such interactions re- 
veal new phenomena in the wave dynamics in discrete 
systems, and due to their potential applications to all- 
optical switching (see Ref. [1] for a review of the light 
transmission in linear and nonlinear discrete systems). 
It was predicted that such lattices support discrete soli- 
tons [2, 3, 4], i.e., localized modes for which the discrete 
diffraction in the waveguide array balances the onsite cu- 
bic (x^) nonlinearity, and more complex modes, such as 
discrete vortex solitons [5] (in two-dimensional lattices), 
discrete surface solitons [6] and light bullets [7] . In addi- 
tion, the existence of continuous solitons in semi-discrete 
nonlinear media based on Bragg optical fibers have also 
been recently predicted [8]. The discrete solitons, as well 
as their counterparts predicted in x^ media [9], have 
been observed in experiments [10, 11, 12]. 

In recent work [13], it was demonstrated that, if a peri- 
odic array of optical waveguides is nonlinearly coupled to 
a slab waveguide, both structures being made of a mate- 
rial with quadratic (x ) nonlinearity, the combined op- 
tical structure supports a new species of optical solitons, 
namely semi-discrete composite solitons (SDCS). These 
solitons are different from their fully discrete and con- 
tinuous counterparts in several aspects; in particular, it 
has been shown that both even (symmetric) and twisted 
(antisymmetric) SDCS may be stable. 

In this paper, we extend those ideas to the obviously 
relevant case of composite structures made of optical me- 
dia with cubic nonlinearity. We demonstrate, for the first 
time to our knowledge, that such media support onc- 
and two-parameter families of stable SDCSs, which con- 
sist of two mutually trapped components, viz., a discrete 
one, which is chiefly carried by the waveguide array, and 
a continuous component, bound to the slab waveguide. 
We also demonstrate that such SDCS can be used in opti- 
cally controlled beam steering devices, and therefore they 
could have important practical applications. 

We consider the structure built as a periodic array of 



stripe optical waveguides, whose adjacent cores are sep- 
arated by a distance do, which is coupled to a single- 
mode slab waveguide. The waveguide array may be ei- 
ther buried into the slab waveguide, as shown in Fig. 1, 
or mounted on top of it. Both the stripes and the slab 
are to be made of a Kerr material, possible choices being 
silicon-on-insulator (SOI) [14], polymers [15], or semicon- 
ductor heterostructures [16]. We also assume that the 
structure is designed so as to make the intrinsic modes 
in the stripe and slab waveguides polarized orthogonally 
to each other, thus ruling out the linear coupling between 
them. Specifically, if we consider buried SOI waveguides, 
by simply changing the transverse dimensions of the sil- 
icon waveguide one can design stripe waveguides which 
support either TE or TM polarized modes [17], and there- 
fore an experimental set up involving orthogonally po- 
larized modes can be readily implemented. In the ex- 
periment, the incident wide beam may be unpolarized, 
the stripe and slab waveguides picking mutually orthog- 
onal polarization components from it. As a result, the 
propagating mode in each stripe is linearly coupled to its 
counterparts in the adjacent ones, and nonlinearly cou- 
pled, through the cross-phase modulation (XPM), to the 
slab mode. In addition, the case of linear coupling be- 
tween a transversely confined waveguide mode and a slab 
waveguide mode has previously been considered in Ref. 
[18] and therefore it would not be discussed here. 

Coupled-mode equations for optical fields in the com- 
posite medium can be derived in a consistent form 
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FIG. 1: (Color online) A schematic of the composite optical 
structure. 



2 



[19, 20], resulting in a system of ordinary differential 
equations for the discrete component, coupled to a par- 
tial differential equation governing the propagation of the 
slab mode, 
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Here, (3 Ctd = (3 c ,d/cd are normalized wavenumbers, with 
0c (Pd) the propagation constant in the slab (stripe) 
waveguide in physical units, and Cd the constant of 
the linear coupling between the stripe modes, while 
£ = z/zq and r\ — x/xq , where z and x are the lon- 
gitudinal and transverse coordinates, and zq = l/cd, 
xq = (/? c Cd)~ 1//2 . The normalized fields 4> n and * are 
rescaled modal fields u n and U (which are measured in 
units of y/W and ^/W/m, respectively) : 4> n — \/"fd/cdU n 
and \V = \J"ic/ c dU ■ Finally, Hd and k c are relative 
strengths of the XPM and SPM (self-phase modulation) 

interactions, n d = y/jdc/lc k c = \/^ c d/ld- Here, co- 
efficients 7 are defined as integrals over cross sections 
of the respective waveguides (Ad and A c , the latter de- 
fined as the transverse area of the slab waveguide be- 
tween adjacent stripes), involving the convolution of the 
third-order susceptibility tensor, x*- 3 -*, with the stripe 
and slab modes, ed(x,y) and e c (y), which are normal- 
ized to the corresponding mode powers, Vd and V c (mea- 
sured in W and W/m): 7d = (3e u;/4Vj) J A dAe* d ■ 
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(3e u/2VdV c d )f Ats dAe c • x t3, :e c e^e rf . The same system 
of coupled-mode equations may also be interpreted as the 
one governing the transmission of waves <p n and \V with 
different carrier frequencies, which also rules out the lin- 
ear coupling between them. 

To focus on the most fundamental case, we assume 
that the stripe and slab modes are phase-matched, (3d = 
f3 c = (3, hence the corresponding linear terms in Eqs. (1) 
can by removed (this condition can be easily satisfied by 
adjusting geometrical parameters of the stripe and slab), 
and the nonlinear coupling between the stripe and slab 
waveguides is symmetric, with K d = Kc = ^- 

Then, Eqs. 

(1) conserve the Hamiltonian, 
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FIG. 2: (Color online) Field profiles of the three species of the 
semi-discrete composite solitons: odd, alias onsite-centered 
(a); even, alias intersite-centered (b); and twisted (c). 
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FIG. 3: (Color online) Top panels: total power P and power 
Pd of the discrete component (solid and dashed lines, respec- 
tively) versus A, for (a) odd and (b) even solitons. (c) Hamil- 
tonian H of the odd and even solitons (solid and dashed lines, 
respectively) versus Pd- (d) The dominant instability growth 
rate for even solitons versus A. 



(the asterisk stands for the complex conjugation), and 
total powers of the discrete and continuous components, 
p d = X)„ l^nl 2 and Pc = l^fdr], respectively. 

General soliton solutions to Eqs. (1) can be looked 
for as <f> n (0 = <f> n e iXdC and ^(f],() = ^(il)e iXaC , with 
Ad 7^ A c , i.e., they form a two-parameter family of soli- 
tons [21, 22, 23]. To find these solutions, we numeri- 
cally solved the respective stationary version of Eqs. (1), 
thus constructing both a particular one-parameter set of 
the solutions, with Ad = A c = A, and the general two- 
parameter family of SDCSs (Ad 7^ A c ). 

Figures 2(a)-(c) display generic examples of the soli- 
tons belonging to the one-parameter families of the odd, 
even, and twisted types, obtained through this procedure 
(the former two may also be classified as onsite-centered 
and intersite-centered, respectively, as concerns their dis- 
crete component). The entire families arc presented in 
Figs. 3(a)-(c) by plots showing Hamiltonian H, total 
power P = P c + Pd, and the power in the discrete com- 
ponent, Pd, as functions of common wavenumber A. Both 
odd and even solitons exist only if A exceeds a threshold 
value, Athr, which increases with the coupling strength, 
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FIG. 4: (Color online) Unstable propagation of even (top pan- 
els) and twisted (bottom panels) semi-discrete composite soli- 
tons. Left and right panels correspond to the discrete and 
continuous components, respectively. The soliton parameters 
are the same as in Fig. 2. 



k. This result can be explained by the dispersion rela- 
tion for CW solutions to Eq. (lb), with amplitude W 
and transverse wavenumber ko, A cw = 2 cos fco + 
since nonlinear modes exist for A > A cw only, A t hr indeed 
increases with k. Moreover, diagram (H, Pd) in Fig. 3(c) 
suggests that the odd (onsite-centered) composite soli- 
tons may be stable, as they correspond to a smaller value 
of H, and become more localized as Pd increases. Indeed, 
note that the Peierls-Nabarro barrier (PNB), which is de- 
fined as the difference in H of odd and even solitons that 
have the same power Pd [24], increases with Pd. 

The expectations concerning the stability of the SD- 
CSs are supported by the linear-stability analysis. To 
perform the analysis, we linearized Eqs. (1) around 
the stationary soliton solutions, 4 > ^°^ and vi/( so1 ), set- 
ting 0„ = $f ol) + (<50„e J ^ + 5j>* n e- iv *<)e iX t and * = 
#(sol) + (s^e^C + S^*e- iu *< : )e tX(: , and calculating eigen- 
values v from the ensuing linearized equations for the 
small perturbations. The results, summarized in Fig. 
3(d), demonstrate that all even solitons are unstable, 
which is natural to states whose discrete component is 
centered at the intersite position [1]. Twisted SDCSs are 
unstable too, in their entire domain of existence. On the 
contrary, odd SDCSs are completely stable, as Im(f) = 
for all A at which they exist. 

A generic scenario for the unstable propagation of even 
and twisted SDCSs is illustrated in Fig. 4, the input 
solitons being those presented in Fig. 2. Thus, in the case 
of even SDCSs, the discrete component evolves into an 
odd discrete soliton, which then becomes locked in with 
the continuous component. As a result, an odd SDCS is 
formed. Moreover, during this process part of the energy 
contained in the initial soliton is released as radiation. 
On the other hand, twisted SDCSs break in a pair of 
mutually trapped odd solitons, which continuously emit 
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FIG. 5: (Color online) Field profiles of the three species of 
the semi-discrete composite solitons: odd, (a); even, (b); and 
twisted (c), corresponding to the case k < 0. 




FIG. 6: (Color online) Unstable propagation of odd (top pan- 
els), even (middle panels), and twisted (bottom panels) semi- 
discrete composite solitons. Left and right panels correspond 
to the discrete and continuous components, respectively. The 
soliton parameters are the same as in Fig. 5. 



radiation as they copropagate. 

For the sake of completeness of the study of the model 
(1), we also investigated the existence and stability of 
soliton solutions in a more formal case, namely that of 
k < 0. Generic odd, even, and twisted soliton solutions 
are shown in Fig. 5. Our analysis shows that the main 
difference between this case and the case in which k > 
is that for k < all three species of solitons are un- 
stable. Thus, as our numerical simulations illustrate, if 
k < odd SDCSs are unstable, too, upon a small trans- 
verse perturbation in the field profile. Specifically, due 
to the self-defocusing characteristics of the XPM interac- 
tion in this case, the discrete and continuous components 
repel each other and form discrete and continuous soli- 
tons, which separate from each other upon propagation 
(see Fig. 6, top panels). This soliton break up is also 
observed in the case of even and twisted SDCSs. 

General two-parameter SDCS families of the odd and 
even types, with A^ ^ A c , have also been found numer- 
ically as stationary localized solutions of Eqs. (1). The 
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FIG. 7: (Color online) The power in the discrete and con- 
tinuous components (middle and top panels, respectively) of 
the general soliton family versus independent wavenumbers 
Ad and A c . The dependence of Hamiltonian H on Ad and A c 
is shown in the bottom panel. In all panels, the nonlinear- 
coupling constant is k = 0.1, the blue and green (upper and 
lower) surfaces corresponding to the odd and even SDCS, re- 
spectively. 

results of these calculations, performed for both odd and 
even SDCSs, are presented in Fig. 7, through the respec- 
tive dependences of the dynamical invariants, P dl P c and 
H, on X d and A c . The first among these dependences 
shows that, similar to the case of A^ = A c = A, SDCSs 
exist only if X d exceeds a certain threshold value. 

To analyze the stability of the general SDCS family, we 
first used an extension of the Vakhitov-Kolokolov (VK) 
criterion [25] for two-parameter families of solitons [26, 
27, 28], according to which the stability changes across 
curves J = in the parameter space, where the Jacobian 
is given by the following expression: 

= d(P d , P c ) = dP d dP c dP d dP c 

d(\ d ,x c ) d\ d d\ c ox c dx d u 

Using the computed powers P d and P c , we have con- 
cluded that J 7^ (in fact, J > 0) for all values of A^ and 
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FIG. 8: (Color online) Contour plots of the discrete and con- 
tinuous field components (left and right panels, respectively) 
of odd (stable) solitons subjected to the initial phase tilt, fco- 
In this and following figures, the profiles were recorded after 
propagation distance £o = 10 , for coupling constant n — 0.1. 
At Pd < 1, the discrete component is spread over several 
stripe waveguides, a transition to pinning being clearly seen 
at Pd ~ 2 and 3, for fco = 0.4 and 0.7, respectively. Note 
that the continuous component is stronger localized than the 
discrete one. 



A c at which SDCSs exist; therefore, both odd and even 
SDCSs do not change the VK stability in their existence 
domains. 

Because the VK criterion is only a necessary stability 
condition [note that it does not predict the instability of 
the one-parameter family of even SDCSs, see Fig. 3(b)], 
we have also performed the full linear-stability analysis 
for the two-parameter families, and concluded that (as 
might be expected), the odd solitons are stable, whereas 
their even counterparts are not. These conclusions ac- 
cord to the topology of surface H(P d , P c ), as seen in Fig. 
7(c): since, for both the odd and even SDCS families, 
H(P d , P c ) is a smooth single- valued surface with no cus- 
pidal edges (Whitney surfaces), neither type of the soli- 
tons should change its stability within its existence do- 
main [27, 29]. 

The enhanced localization of the composite solitons 
with increasing power P d of the discrete component can 
be exploited in all-optical beam-steering devices [4, 9, 30]. 
To demonstrate this possibility, we consider SDCSs with 
a phase tilt , accounted for by transverse wavenumber fco > 
at the input facet of the waveguide array, {4> n , ^} lc=o — 

|^n° l \ 1 jK so1 ) J e ife °'', and follow their propagation in the 

waveguide array. Figure 8 presents the results for odd 
(stable) solitons, in the case of X d = X c = X. As might 
be expected, at low power P d the "kicked" (phase-tilted) 
SDCS readily moves across the waveguide array. How- 
ever, as P d increases, the transverse displacements of 
both the discrete and continuous components decrease, 
and, for P d exceeding a certain critical value [in Fig. 8, 
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FIG. 9: (Color online) The location of the peak intensity in 
the discrete and continuous components (left and right panels, 
respectively) of the odd soliton versus power Pd of its discrete 
component, at different values of phase pitch fco applied at 
C = 0. The coupling constant is k = 0.1. 

(Pd) cr ~ 2 and 3, for k = 0.4 and 0.7, respectively], 
the soliton remains trapped at its initial location. This 
finding is consistent with the fact that the PNB increases 
with Pd, as per Fig. 3(c). 

Dependencies of the transverse displacement of the odd 
SDCS on phase pitch fco, which causes its displacement, 
and power Pd of its discrete component, are displayed in 
Fig. 9. Naturally, the transverse displacement increases 
with fco, as well as critical power (Pd) CI required to pin 
the soliton. It is worthy to note that the discrete and 
continuous components remain mutually trapped even 
close to the linear limit (at small Pd), which is not obvious 
in Fig. 8. 

To summarize, we have demonstrated that the struc- 



ture built as a single-mode slab waveguide coupled to 
a discrete waveguide array, with both parts made of 
a Kerr-nonlinear material, supports semi-discrete com- 
posite solitons of odd (onsite-centered) , even (intersite- 
centered), and twisted types. One- and two-parameter 
families of the odd solitons (with one or two independent 
wavenumbers, respectively) are stable, while the other 
types are not. Note that in the case of optical media with 
quadratic nonlinearity both even and twisted SDCSs may 
be stable. This should not be a surprise as quadratically 
and cubically nonlinear media represent essentially dif- 
ferent optical system and therefore one expects that they 
support nonlinear modes with different physical proper- 
ties. In addition, in the case of quadratically nonlinear 
media, SDCSs have a more complex structure, as the dis- 
crete and the continuous components can be excited at 
either the fundamental frequency or at the second har- 
monic. 

The potential application of the stable solitons, kicked 
in the transverse direction, to all-optical beam steer- 
ing was also demonstrated. The analysis presented here 
for the one-dimensional setting can be extended to two- 
dimensional lattices and photonic-crystal fibers. For ex- 
ample, nonlinear modes supported by 2D nonlinear opti- 
cal media invariant to discrete symmetry transformations 
have already been investigated [31, 32, 33, 34]. It may 
also be interesting to analyze semi-discrete solitons in 
media with the self-defocusing cubic nonlinearity, where 
the discrete component would be of the staggered type. 

This work was supported by NSF under Grant No. 
ECS-0523386, and by DARPA/AFOSR under Grant No. 
FA9550-05- 1-0428. 
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